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D , Abstract 

The elastic electroproduction of p mesons is studied at HERA with the HI detector for 
a photon virtuality in the range 1 < Q 2 < 60 GeV 2 and for a hadronic centre of mass 
energy in the range 30 < W < 140 GeV. The shape of the (tttt) mass distribution in 
the p resonance region is measured as a function of Q 2 . The full set of p spin density 
matrix elements is determined, and evidence is found for a helicity flip amplitude at the 
level of 8 ± 3% of the non-flip amplitudes. Measurements are presented of the dependence 
of the cross section on Q 2 , W and t (the four-momentum transfer squared to the proton). 
They suggest that, especially at large Q 2 , the 7*p cross section develops a stronger W 
dependence than that expected from the behaviour of elastic and total hadron— hadron cross 
sections. 
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1 Introduction 

Measurements of the elastic electroproduction of vector mesons at HERA over a wide range of 
exchanged photon virtuality Q 2 are of particular interest. For many years it has been known that 
at low Q 2 , that is with no hard scale, vector meson electroproduction exhibits all the properties 
of a soft diffractive process. Predictions of soft processes based on QCD calculations are how- 
ever intractable. The presence of a hard scale, that is a significant Q 2 , makes perturbative QCD 
calculations possible. Measurements of the Q 2 dependences of observables in vector meson 
electroproduction thereby provide insight into the transition and the interplay between soft and 
hard processes in QCD. 

This paper presents an analysis of elastic p meson electroproduction: 

e + p — > e + p + p , p-^7r + + 7r~, (1) 

in the Q 2 range from 1 to 60 GeV 2 (Q 2 = — q 2 , where q is the four-momentum of the in- 
termediate photon) and the W range from 30 to 140 GeV (W is the hadronic centre of mass 
energy). 

The data were obtained with the HI detector in two running periods of the HERA collider, 
operated with 820 GeV protons and 27.5 GeV positrons.|] A low Q 2 data set (1 < Q 2 < 5 GeV 2 ) 
was obtained from a special run in 1995, with the ep interaction vertex shifted by 70 cm in 
the outgoing proton beam direction; it corresponds to an integrated luminosity of 125 nb -1 . 
A larger sample with 2.5 < Q 2 < 60 GeV 2 was obtained in 1996 under normal running 
conditions; it corresponds to a luminosity of 3.87 pb _1 . 

The present measurements provide detailed information in the region 1 < Q 2 < 8 GeV 2 and 
they increase the precision of the HI measurement of p electroproduction with Q 2 > 8 GeV 2 , 
which was first performed using data collected in 1994 [[[p. They are compared to results of the 
ZEUS experiment [J2|] at HERA and of fixed target experiments [Q-§]] . 

The HI detector, the definition of the kinematic variables and the event selection are intro- 
duced in section ^. Acceptances, efficiencies and background contributions are discussed in 
section [3[ The shape of the (nn) mass distribution and the evolution with Q 2 of the skewing of 
this distribution are studied in section |j. Section |5| is devoted to the study of the p meson decay 
angular distributions and to the measurement of the 15 elements of the spin density matrix, as 
a function of several kinematic variables. The Q 2 dependence of the ratio R of the longitudinal 
to transverse YP cross sections is measured. The violation of s-channel helicity conservation, 
found to be small but significant at lower energies [0,@], is quantified. Finally, section || presents 
the t distribution and the measurement of the 7*p — » pp cross section as a function of Q 2 and 
W. Predictions of several models are compared to the measurements in sections |] and |[ 



2 HI Detector, Kinematics and Event Selection 

Events corresponding to reaction ([[]) are selected by requiring the detection of the scattered 
electron and of a pair of oppositely charged particles originating from a common vertex. The 

1 In the rest of this paper, the word "electron" is generically used for electrons and positrons. 
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absence of additional activity in the detector is required, since the scattered proton generally 
escapes undetected into the beam pipe. 

HI uses a right-handed coordinate system with the z axis taken along the beam direction, 
the +z or "forward" direction being that of the outgoing proton beam. The x axis points towards 
the centre of the HERA ring. 



2.1 The HI Detector 

A detailed description of the H 1 detector can be found in ^ . Here only the detector components 
relevant for the present analysis are described. 

The scattered electron is detected in the SPACAL [j|], a lead - scintillating fibre calorimeter 
situated in the backward region of the HI detector, 152 cm from the nominal interaction point. 
The calorimeter is divided into an electromagnetic and a hadronic part. The electromagnetic 
section of the SPACAL, which covers the angular range 153° < 9 < 177.5° (defined with 
respect to the nominal interaction point), is segmented into cells of 4 x 4 cm 2 transverse size.0 
The hadronic section is used here to prevent hadrons from being misidentified as the scattered 
electron. In front of the SPACAL, a set of drift chambers, the BDC, allows the reconstruction 
of electron track segments, providing a resolution in the transverse direction of 0.5 mm. 

The pion candidates are detected and their momentum is measured in the central tracking 
detector. The major components of this detector are two 2 m long coaxial cylindrical drift 
chambers, the CJC chambers, with wires parallel to the beam direction. The inner and outer 
radii of the chambers are 203 and 451 mm, and 530 and 844 mm, respectively. In the forward 
region, the CJC chambers are supplemented by a set of drift chambers with wires perpendicular 
to the beam direction. The measurement of charged particle transverse momenta is performed in 
a magnetic field of 1 . 1 5 T, uniform over the full tracker volume, generated by a superconducting 
solenoidal magnet. For charged particles emitted from the nominal vertex with polar angles 
20° < 9 < 160°, the resolution on the transverse momentum is Ap t /pt — 0.006 p t (GeV). 
Drift chambers with wires perpendicular to the beam direction, situated inside the inner CJC 
and between the two CJC chambers, provide a measurement of z coordinates with a precision 
of 350 fim. 

The (x, y, z) position of the interaction vertex is reconstructed for each event by a global 
fit of all measured charged particle trajectories. For each electron fill in the accelerator, a fit is 
performed of the dependence on z of the mean x and y positions of the vertices. This provides 
a measurement of the corresponding beam direction, which varies slightly from fill to fill. 

The absence of activity in the HI detector not associated with the scattered electron or the p 
decay is checked using several components of the detector. The liquid argon (LAr) calorimeter, 
surrounding the tracking detector and situated inside the solenoidal magnet, covers the polar 
angular range 4° < 9 < 154° with full azimuthal acceptance. The muon spectrometer (FMD), 
designed to identify and measure the momentum of muons emitted in the forward direction, 
contains six active layers, each made of a pair of planes of drift cells, covering the polar angular 
region 3° < 9 < 17°. The three layers situated between the main calorimeter and the toroidal 

2 In this paper, "transverse" directions are relative to the beam direction. 



magnet of the FMD can be reached by secondary particles arising from the interaction of small 
angle primary particles hitting the beam collimators or the beam pipe walls. Secondary particles 
or the scattered proton at high \t\ can reach a set of scintillators, the proton remnant tagger 
(PRT), placed 24 m downstream of the interaction point and covering the angles 0.06° < 9 < 
0.17°. 



2.2 Kinematic Variables 

The reconstruction method for the kinematic variables has been optimised for the p measure- 
ment. 

The Q 2 variable is computed from E , the incident electron beam energy, and the polar 
angles 9 e and 9 P of the electron and of the p meson candidates [§]: 

AE 2 

O 2 = --2 (2) 

* tan(0 e /2) (tan(0 e /2) + tan(0„/2) ) ' V 

The electron emission angles are determined using the reconstructed vertex position and the 
track segment in the BDC corresponding to the electron cluster candidate. The momentum of 
the p meson is reconstructed as the sum of the momenta of the two pion candidates: 

Pp = Pn+ + V-k- ■ (3) 

The inelasticity y is defined as 

V = r , (4) 

p ■ k 

where p and k are the four-momenta of the incident proton and of the incident electron, respec- 
tively. For this analysis, y is computed, with very good precision, using the energy, E p , and the 
longitudinal momentum, p Zp , of the p meson candidate (JT0[] : 

E — Pz 

v = ^ieV (5) 

The hadronic mass, W, is computed using the relation 

W 2 = ys - Q 2 , (6) 

where s is the square of the ep centre of mass energy. 

The variable t is the square of the four- momentum transfer to the proton. At HERA energies, 
to very good precision, its absolute value is equal to the square of the transverse momentum 
of the outgoing proton. The latter is computed, under the assumption that the selected event 
corresponds to reaction ([]]), as the sum of the transverse momenta p t of the p meson candidate 
and p te of the scattered electron: 

1 2 



\Pt P +PtX- (7) 



The value of t is thus distorted if the event is due to the production of a hadron system of which 
the p is only part and of which the remaining particles were not detected. For use in eq. (0), p te 
is determined from the p candidate measurement and the electron beam energy, such that 

_ 2E -E p + p Zp 
Vu ~ tan(0 e /2) • (8) 

This relation assumes reaction (|T]) and the absence of QED radiation. 

Finally, the total event E — p z variable is computed as the sum of the differences between the 
energies and the longitudinal momenta of the electron and pion candidates, where the electron 
energy measured in the SPACAL calorimeter is used. 



2.3 Trigger and Event Selection 

The trigger and selection criteria for the events used in this analysis are summarised in Table [I]. 
Events are selected only from runs for which all relevant parts of the detector were functioning 
efficiently. 

For the 1995 shifted vertex run, the trigger was based on the detection of a cluster in the 
electromagnetic section of the SPACAL calorimeter with energy greater than 12 GeV. For the 
1996 data, the energy threshold was increased to 15 GeV and, in order to reduce the rate of 
background events due to synchrotron radiation from the electron beam, the centre of gravity of 
the cluster was required to lie outside the innermost part of the SPACAL, with — 16 < x < 8 cm 
and — 8 < y < 16 cm. Independent triggers were used to determine the efficiency of this trigger. 

Off-line, electron candidates are defined as well identified electromagnetic clusters in the 
SPACAL with energy larger than 17 GeV, correlated with a track segment in the BDC. The 
transverse position of the BDC track segment has to be more than 8.7 cm from the beams for 
the 1995 data sample, and must correspond to the region of the SPACAL included in the trigger 
for the 1996 data. 

Exactly two oppositely charged pion candidates are required, with polar angles of emission 
20° < 9 < 160° (1996 data []), and transverse momenta with respect to the beam direction 
p t > 0.1 GeV, so that detection and reconstruction in the central tracker are efficient. The 
reconstructed interaction vertex has to lie within 30 cm in z of the nominal interaction point. 

Rejection of p meson events with proton dissociation and of other backgrounds is achieved 
using three selection criteria: there must be no cluster in the LAr calorimeter with energy greater 
than 0.5 GeV that is not associated with the pion candidates, there must be no more than one hit 
pair recorded in the FMD and there must be no signal in the PRT Given the limiting angle of 
20° for pion candidates, this corresponds to requiring no activity for a range in pseudorapidity 
1.75<7?<7.5f3 

The cuts 1 < Q 2 < 60 GeV 2 and 30 < W < 140 GeV, which define the kinematic domain 
under study, correspond to the region in which the electron and hadronic track acceptances are 

3 For the 1995 data, no cut on the track polar angle is made. 

4 The pseudorapidity rj of an object detected with polar angle 6 is defined as rj = — In tan(0/2). 



Trigger 


cluster in SPACAL with energy > 12 (15) GeV in 1995 (1996) 
and with -16 < x < 8 cm and -8 < y < 16 cm (1996) 


Electron 


cluster in electromagnetic SPACAL with energy > 17 GeV 
distance between cluster c.o.g. and BDC track < 3 cm 
BDC segment > 8.7 cm from the beams (1995) 
transverse width of cluster < 3.2 cm 
energy in hadronic SPACAL < 0.2 GeV 


Pion candidates 


exactly two tracks with opposite signs 

20° < 6 < 160° (1996) 

particle transverse momenta p t > 0.1 GeV 

vertex reconstructed within 30 cm of nominal position in z 


Additional activity 


no cluster in LAr with energy > 0.5 GeV 
at most 1 hit pair in FMD 
no hit in PRT 


Mass selection 


0.6 < m™ < 1.1 GeV 
rriKK > 1-040 GeV 


Kinematic domain 

1995 data 

1996 data 


1.0 < Q 2 < 5 GeV 2 , 40 < W < 140 GeV 
2.5 < Q 2 < 4 GeV 2 , 30 < W < 100 GeV 
4.0 < Q 2 < 6 GeV 2 , 30 < W < 120 GeV 
6.0 < Q 2 < 14 GeV 2 , 40 < W < 140 GeV 
14.0 < Q 2 < 60 GeV 2 , 50 < W < 140 GeV 


Other cuts 


|£| < 0.5 GeV" 
E-p z >Ah GeV 



Table 1 : Summary of trigger conditions and event selection criteria (see text for details). 



high. A cut |i| < 0.5 GeV 2 is also applied, the purpose of which is threefold. Firstly, the 
acceptance for elastic events decreases at larger \t\ values, because the probability becomes sig- 
nificant that the proton hits the beam pipe walls, thus producing a signal in the PRT. Secondly, 
the \t\ cut suppresses events from processes which are not elastic and have a flatter t distribu- 
tion, in particular p production with proton dissociation. Thirdly, it suppresses the production 
of hadron systems of which the p is only part and in which the remaining particles were not 
detected, thereby distorting the measurement of t (see eq. [7]). A further cut, E — p z > 45 GeV, 
is designed to minimise the effects of initial state photon radiation from the electron. 

The selected domain for m n77 , the invariant mass of the two pion candidates, is restricted 
to 0.6 < m^ < 1.1 GeV, which covers the p meson mass peak and avoids regions with large 
background contributions. In order to minimise 4> meson contamination, the invariant mass 
of the pion candidates is also computed with the assumption that they are kaons, and the cut 
m KK > 1.040 GeV is applied on the corresponding m KK mass. 

After all selection cuts, the 1995 sample (1 < Q 2 < 5 GeV 2 ) contains about 500 events, 
and the 1996 sample (2.5 < Q 2 < 60 GeV 2 ) 1800 events. 



3 Detector Effects and Background Contributions 
3.1 Acceptances and Efficiencies 

Acceptances, efficiencies and detector resolution effects are determined using the DIFFVM 
Monte Carlo simulation [JTTJ] , a program based on Regge theory and the vector meson domi- 
nance model (VDM). The simulation parameters are adjusted following the measurements pre- 
sented below for the dependence of the cross section on Q 2 , W, t and for the p meson angular 
decay distributions. The detector geometry and its response to generated particles are simulated 
in detail. The same reconstruction procedures and event selection criteria as for real events are 
applied. As an illustration of the good quality of the simulation, Fig. [I] presents a comparison 
of the distributions of several variables for the data and for the Monte Carlo simulation. The 
distribution of the azimuthal angle of the p meson (Fig. [j]c) reflects the regions of the SPACAL 
that are active in the trigger. The distribution of the transverse momenta of the pion candidates 
(Fig. [TJd) depends on the details of the p meson decay angular distribution. It has been carefully 
checked that the Monte Carlo simulation reproduces well the details of the tracker acceptance 
and efficiency, both for positively and for negatively charged pions. 

In the kinematic domain defined in Table [j], the acceptance depends most strongly on Q 2 in 
a purely geometrical manner related to the trigger conditions. The cuts on the polar angles and 
on the minimum transverse momenta of the pion candidates induce VT-dependent acceptance 
corrections, which are sensitive to the angular decay distributions. The Q 2 and W limits of the 
selected kinematic domain are such that the efficiency is almost constant over each bin. The cut 
on E — p z induces very small corrections. 

For each of the measurements presented below, systematic errors are computed by varying 
the reconstructed polar angle of the electron by ±0.5 mrad, which corresponds to the systematic 
uncertainty on this measurement, and by varying in the Monte Carlo simulation the cross section 
dependence on Q 2 , W, t and the p meson decay distributions by the amount allowed by the 
present measurements (see Jl~2| ] for more details). Small remaining uncertainties related to the 
simulation of the tracker uniformity are neglected. Further systematic uncertainties that affect 
only certain measurements are described where appropriate below. The positive and negative 
variations are combined separately in the form of quadratic sums, to compute the systematic 
errors. 

In addition to the effects studied with the DIFFVM simulation, the trigger efficiency is stud- 
ied using several independent triggers. Regions of the SPACAL for which the trigger efficiency 
is below 94% are discarded from the measurement. Losses of elastic events due to noise in 
the LAr, FMD and PRT detectors are computed from randomly triggered events in the detector. 
Radiative corrections are determined using the HERACLES program [|13[j. 



3.2 Background Contributions 

The main background contributions to p meson elastic production are due to the elastic produc- 
tion of uj and (ft mesons and to diffractive p production with proton dissociation. 
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Figure 1 : Uncorrected distributions of the polar angle of the scattered electron, the polar angle 
of the p meson, the azimuthal angle of the p meson in the laboratory frame, and the transverse 
momenta of the two pion candidates, for the 1996 data sample (points) and for the Monte Carlo 
simulation (histograms), after all selection cuts. 
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3.2.1 Elastic Production of uj and </> Mesons 

The elastic production of a; mesons: 

e + p — > e + uj + p (9) 

may produce background in the present data sample through the two uj decay modes [IT3] : 



uj -»• 7T+ + 7T" + 7T° (BR = 88.8%) (10) 

W ->• 7T+ + 7T _ (BR = 2.2%) . (1 1) 

The contribution of the first decay mode is efficiently reduced by the mass selection cut, by 
requiring the absence in the LAr calorimeter of clusters with energy larger than 0.5 GeV which 
are not associated with a track, and by the cut on the variable t. However, events from the second 
decay mode are selected within the present sample. This background is subtracted statistically 
assuming the uj : p ratio of 1 : 9 which is motivated by SU(3) flavour symmetry and is consistent 
with HERA photoproduction measurements |JT3|] . 

The production rate of mesons: 

e + p — > e + <fi + p (12) 

amounts to about 15% of the p production rate for the present kinematic domain [TO-ffBfl. The 



following decay modes [ [14| ] may lead to the presence of background events in the selected 
sample: 

cf)^K + + K- (BR = 49.1%) (13) 

</> -»• p + it (BR = 12.9%) (14) 

<f) ->• 7T+ + 7T- + 7T° (BR = 2.7%) . (15) 

The first contribution is mostly eliminated by the m KK and the m^ mass selection cuts, and 
the other two are significantly reduced by the cuts against additional particles and by the t and 
mass selection cuts. 

Using the DIFFVM Monte Carlo simulation, the contribution of uj and <p elastic production 
remaining in the selected sample is determined to be 3.3 ± 2.0% in the invariant mass range 
0.6 < trimr < 1.1 GeV, where 1.4% and 1.9% come from the uj and contributions, respec- 
tively. For the study of the shape of the mass distribution, the m^ range used is 0.5 < m niT < 
1.1 GeV, where the contributions of a; and 4> elastic production are determined to be 4.7% and 
2.3%, respectively, and are subtracted statistically bin-by-bin from the mass distributions (see 
section HI). 



3.2.2 Diffractive Production of p Mesons with Proton Dissociation 

An important background to elastic p production is due to the diffractive production of p mesons 
with proton dissociation 

e + p^ e + p + Y (16) 
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when the baryonic system Y is of relatively low mass M Y < 1.6 GeV and its decay products 
are thus not detected in the PRT, the FMD or the forward regions of the LAr calorimeter and 
the tracking detector. 

The contamination from proton dissociation is determined using the DIFFVM Monte Carlo. 
The distribution of M Y is generated as (see [|T9|]): 



(IMy My ' 



For M Y < 1.9 GeV, the details of baryonic resonance production and decays are simulated fol- 
lowing the Particle Data Group (PDG) tables Jl4t|. For larger masses, the system Y is modelled 
as formed of a quark and a diquark, which fragment according to the JETSET algorithm pO]. 



The t distribution of proton dissociation events is modelled by an exponentially falling distri- 
bution with a slope parameter b = 2.5 GeV -2 (cf. the measurements in [ |T7| ] and [|ZT|]). The 
DIFFVM Monte Carlo is also used to compute the probability that the scattered proton in an 
elastic p event with \t\ < 0.5 GeV 2 gives a signal in the PRT. 

The proton dissociation background in the selected sample of events is determined without 
making any hypothesis for the relative production rates for elastic and inelastic events. It is 
deduced using the total number of p events and the number of p events with no signal in the 
PRT or the FMD, given the probabilities of obtaining no signal in these detectors for elastic 
interactions and for interactions with proton dissociation. These probabilities are determined 
using the Monte Carlo simulation. The proton dissociation background in the present sample 
amounts to 11 ± 5%. The uncertainty on this number is estimated by varying by ±0.3 the 
exponent of M Y in eq. (JTVp, by varying the slope parameters of the exponential \t\ distributions 
of elastic and proton dissociation events within the experimental limits (see section |6TT| ) and by 
computing the correction using only the PRT or only the FMD [|12|]. 



3.2.3 Other Background Contributions 

Other background contributions are negligibly small. The background due to the n + n~n°n° 
decay mode of the p'(1450) meson is determined to be only 1 ± 1%, due to the cuts against 
additional particles and the cut on the variable t. The study of the mass distributions presented 
in section |] also indicates that events with photon dissociation into vector mesons other than 
p, ijj and do not contribute more than 1%.Q The background from photoproduction events 
with a hadron being misidentified as the electron candidate in the SPACAL is extremely small, 
because of the high E e cut. 

5 In the analysis of the 1994 data [En, events were accepted with a maximum energy of 1 GeV for clusters in 
the LAr calorimeter which are not associated with tracks. A contribution of 11 ± 6% non-resonant background, 
concentrated mainly at small m^n masses, was thus subtracted from the cross section measurement. For the present 
analysis, the limit on the cluster energy is 0.5 GeV, leading to a small background contribution, but the losses of 



events due to noise in the LAr calorimeter amount to ~ 10%, as estimated using random trigger data (section 3.1 ) 
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4 Mass Distributions 

For the 1996 events passing the selection cuts of Table [j], with (Q 2 ) = 4.8 GeV 2 and (W) = 
75 GeV, the distribution of m n7r , the invariant 7i + n~ mass, is presented in Figs. and [3| for 
five domains in Q 2 . The to and <p background contributions (see section |3.2.1| ) are subtracted 



according to their mass distribution obtained from the DIFFVM Monte Carlo simulation. 

The mass distributions are skewed compared to a relativistic Breit-Wigner profile: enhance- 
ment is observed in the low mass region and suppression in the high mass side. This effect has 
been attributed to an interference between the resonant and the non-resonant production of two 
pions [[22p. In order to extract the contribution of the resonant part of the cross section, two 
different procedures are used. 



Following the phenomenological parameterisation of Ross and Stodolsky Q23|], the m n7T 
distribution is described as: 

— — = f p BW p (m n7T ) (-^) n + f bg , (18) 

where fp is a normalisation constant and 

BW p {m^) - — -j 2 , 2 ; 2~F27 V ^) 



is a relativistic Breit-Wigner function with momentum dependent width p4| ] 



q* 2 



rw = r P (^-) 3 —— — -r= . (20) 



Here, T p is the p resonance width, q* is the pion momentum in the (7r + 7r~) rest frame and q$ is 
this momentum when m n7T = m p . The factor (m p /m n7T ) n in eq. ( |i~8"l ) accounts for the skewing 
of the shape of the p signal. The background term f bg is parameterised using a distribution in 
phase space which includes the effect of the dipion threshold and an exponential fall off: 

fbg = <*i K - 2m^) Q2 e~ a3m - , (21) 

where m„- is the pion mass and a\, a 2 and a 3 are constants. 

With eq. (|T8|), the mass distribution for all selected events with 2.5 < Q 2 < 60 GeV 2 is fitted, 
after subtraction of the to and background contributions, over the range 0.5 < m CT < 1.1 GeV, 
with the parameters f p , m p , T p , n, a±, a 2 and a 3 left free. The resonance mass is found to be 
0.766 ± 0.004 GeV and the width 0.155 ± 0.006 GeV, in agreement with the PDG values of 
0.770 and 0.151 GeV [JT4]]. The fit value of the skewing parameter is n = 1.4 ± 0.2 and the 
background contribution corresponds to 1 ± 1 % of the number of events in the peak. The fit is 
of good quality: xV ndf = 20.3/17. 

For the five Q 2 domains presented in Fig. 0, fits to the form of eq. ( [T8] ) are thus performed 
with the mass and the width of the p meson fixed to the PDG values and assuming the absence of 
non-resonant background (fb g = 0). This leaves two free parameters: the overall normalisation 
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Figure 2: Acceptance corrected m^ mass distributions for the 1996 data sample, after statistical 
subtraction of the remaining uj and </> background contributions, divided into five bins in Q 2 . The 
superimposed curves are the result of fits to skewed relativistic Breit-Wigner distributions using 
the Ross-Stodolsky parameterisation of eq. (|T8|), with the p mass and width fixed at the PDG 
values and assuming no non-resonant background. The solid curves are the results of the fits, 
the dashed curves correspond to the non- skewed Breit-Wigner contributions. The errors on the 
data are statistical only. 
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Figure 3: Same data as in Fig. g, but compared to the Soding parameterisation of eq. (|22|). 
The solid curves are the results of the fits, the dashed curves correspond to the non-skewed 
Breit-Wigner contributions, and the dotted curves correspond to the interferences between the 
resonant and the non-resonant amplitudes. The errors on the data are statistical only. 
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fp and the skewing parameter n. The results of the fits are presented in Fig. 0, the % 2 /ndf values 
being good in all Q 2 bins. 



The data are also analysed using the parameterisation proposed by Soding [[25|], in which 
the skewing of the mass spectrum is explained by the interference of a resonant p — > n + n~ 
amplitude and ap-wave nn Drell-type background term: 

<l A ' "' - = f P BW p (m nn ) + fj I{m^) + f bg , (22) 



dm 



7T7T 

2 ™2 



m . . , _. 
(m; - m^) 2 + mfp T 2 (m^) 

where fj is a constant fixing the relative normalisation of the interference contribution. In 
view of the uncertainty in the phase between the resonant and the non-resonant amplitudes, no 
constraint is imposed on the relative contributions of the background and interference terms. 

The Soding parameterisation also describes well the Q 2 integrated data in the range 0.5 < 
m^ < 1.1 GeV, with values for the resonance mass and width in agreement with the PDG 
values and non-resonant background compatible with zero. For the five selected Q 2 bins, the 
width and the mass of the p meson are thus fixed and fb g is taken to be zero. Fits to the 
normalisation and the skewing parameter fi/f p are again of good quality, and the results are 
presented in Fig. [|. 

Fig. shows the fit values of the skewing parameters as a function of Q 2 , together with the 
results of other measurements in photoproduction [[2~6]-|2~S|] and in electroproduction [[1],@,|5|] . The 



systematic errors are computed as described in section pTT] , and include in addition the effect 
of the variation by 50% of the u and <\> background contributions. The skewing of the mass 
distribution is observed to decrease with Q 2 . No significant W or t dependence of the skewing 
is observed within the data. 
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Figure 4: Q 2 dependence of the skewing parameters for elastic p production: n, for the Ross- 



Stodolsky parameterisation of eq. (|T8|), and fi/fp, for the Soding parameterisation of eq. (|22|). 
For the present measurements (full circles), the inner error bars are statistical, and the full 
error bars include the systematic errors added in quadrature. The other measurements are from 
HI @ and ZEUS [^H in photoproduction, and from HI 10, ZEUS [|] and E665 [@] in 
electroproduction . 
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5 Helicity Study 

5.1 Angular Decay Distributions 

The study of the angular distributions of the production and decay of the p meson gives in- 
formation on the photon and p polarisation states. The decay angles can be defined in several 
reference frames [p9|]. In the helicity system, used for the present measurement, three angles are 
defined as follows (Fig. |5J). The angle 0, defined in the hadronic centre of mass system (cms), is 
the azimuthal angle between the electron scattering plane and the plane containing the p and the 
scattered proton. The p meson decay is described by the polar angle 9 and the azimuthal angle 
(p of the positive pion in the n + n~ rest frame, with the quantisation axis taken as the direction 
opposite to that of the outgoing proton in the hadronic cms. 

The normalised angular decay distribution iy(cos6>, cp, 0) is expressed following the for- 
malism used in pO] ] as a function of 15 spin density matrix elements in the form 

iy(cosfl, <PA) = ^ { ~(1 - C 4 ) + ^(3 C 4 - 1) cos 2 9 

— v2Rer'j ) Q sin 29 cos (p — r^ sin 2 #cos2yj 

— e cos 20 (rju sin 2 6' + rQ cos 2 9 — v2 Re r\ sin 29 cos tp 

— t\_i sin 2 #cos2<^ ) 
— s sin 20 ( v2 Im r 2 sin 29 sin ip + Im rf_ 1 sin 2 9 sin 2<p J 
+ a/2£ (1 + e) cos ( rf x sin 2 9 + r^ cos 2 9 

— v2 Re tj'q sin 26 1 cos <p — rf_ 1 sin 2 ^cos2(y9 j 
+ \/2e (1 + e) sin I v2 Im r\ sin 2^ sin y3 

+ImrJ_ 1 sin 2 sin 2y» J | , (24) 



where e is the polarisation parameter of the virtual photon: 



(25) 



1-2/ + 2/72' 

with (e) w 0.99 in the present data.[] 

The spin density matrix elements correspond to different bilinear combinations of the he- 
licity amplitudes Tx p \ N , t \^\ N for p meson production, where X p and A 7 are the helicities of the 
p and of the photon, respectively, and Xn and Ajv the helicities of the incoming and outgoing 
proton. The upper indices 1 and 2 of the matrix elements refer to the production of p mesons by 
transverse photons, the index 04 corresponds to a combination of transverse and longitudinal 
photons, and the indices 5 and 6 correspond to the interference between p production by trans- 
verse photons and by longitudinal photons. The lower indices of the matrix elements refer to 
the values of the p meson helicity \ p entering the combination of amplitudes. 

6 In general, there are further contributions to the angular decay distribution, which vanish for unpolarised 
leptons and for e = 1 (see [|30|l). 
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Figure 5: Angle definition for the helicity system in elastic p meson production. 



Specific relations between the amplitudes, leading to predictions for the values of several 
matrix elements, follow from additional hypotheses. 

• s-channel helicity conservation 

For the case of s-channel helicity conservation (SCHC), the helicity of the virtual photon 
is retained by the p meson and the helicity of the proton is unchanged: 



' ApA iy/ ,A-y AjV 



rp A A 

A p A JV ;,A 7 Ajv "A P A 7 "Ajy/Ajv 



(26) 



Single and double helicity flip amplitudes then vanish so that (omitting the nucleon helic- 
ities): 



7a p a 7 — Toi — T__ — T _i — T_i — 0, 
T_u = 71-1 = 0, 

and all matrix elements become zero, except five: 



„04 „1 
00' 



r\_ l , Iiari_ 1 , Re r^ , Im rj . 



Furthermore, the following relationships occur between these elements: 



r\_ x = — Im.r1._4, R- e r io = — Im r 



6 

10 • 



(27) 
(28) 

(29) 

(30) 



Natural parity exchange 

Natural parity exchange (NPE) is defined by the following relations between the ampli- 
tudes:!] 



71 



Ap Ajy/ , — A 7 Ajv 



' 1 \A P — A 7 rp 

v J A p A iv ;,A 7 Ajv 



(3D 



7 For unnatural parity exchange, an additional factor ( — 1 ) appears in the right hand side of eq. ( pl| ) p0| ] . 
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In Table 0, the expressions of the matrix elements are given in terms of the helicity am- 
plitudes for two specific sets of assumptions. In column 2, the double helicity flip amplitudes 
Tx_x and T_ n and the single flip amplitudes T 10 and T_ 10 for the production of transversely 
polarised p mesons by longitudinal photons are neglected, and the NPE relations T _i = — T m 
and T_i_i = T\\ are assumed (see the discussion in section p3J and the presentation of the 



QCD model [pip, in particular eq. d45|), in section |5.4.4| ). In column 3, the matrix elements are 
given for the case of SCHC (i.e. neglecting all helicity flip amplitudes) and assuming the NPE 
relation T_i_i = T n . The nucleon helicities \ N and X N i are omitted from the amplitudes, T, 
for brevity. 

The matrix elements can be measured as the projections of the decay angular distribution 
(eq. ^) onto orthogonal trigonometric functions of the angles 9, ip and </>, which are listed in 
Appendix C of ref. [ pOft . The average values of these functions, for the 1996 data and for the 
kinematic domain defined in Table [U provide the measurements presented in Table ^. The 
results are also presented in Figs. |6|— [8| (and in Tables |j— §) as a function of Q 2 , W and t. 
Statistical and systematic errors are given separately, the systematic errors being computed 
here, and in the rest of section |5j, as described in section pTT] . The data sample is not corrected 
for the small backgrounds due to proton dissociation,^ u and <\> production and radiative effects. 

Within the measurement precision, the matrix elements presented in Table ^] and in Figs. |5J— 18| 
generally follow the SCHC predictions (with the NPE relation T_i_i = T n ). This is not the 
case, however, for the Tq element, which is significantly different from zero (see also the dis- 
cussion of the distribution of the angle 4> in section 153b. It has been checked that this effect is 



not an artifact of the Monte Carlo simulation used to correct the data for detector acceptance 
effects |]I2p. 



As will be discussed in section 1531 the violation of SCHC is small. Information on the 



photon polarisation can thus be obtained from the measurements of the spin density matrix el- 
ements using SCHC as a first order approximation. This analysis is performed in section $72\ 
The violation of SCHC is then studied in more detail in section |53J. For these analyses, the 
good description of the data provided by the function W(cos9, if, <p) is verified through vari- 
ous angular distributions. Finally, section |5J4| presents comparisons of the results with model 
predictions. 



The measurements in [ 17 ] and | T2l| ] indicate that, within errors, elastic and proton dissociation events have the 



same p meson decay angular distributions. 
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Table 2: Spin density matrix elements for the elastic electroproduction of p mesons, expressed 
as a function of the helicity amplitudes T\ p a 7 : second column: the single-flip T 10 and double- 
flip Ti_! amplitudes are neglected and the NPE relations ( CT ) are assumed for the other ampli- 
tudes; third column: the SCHC conditions and the NPE relation T!_i_i = Xii are assumed (i.e. 
the T i helicity flip amplitude is also neglected). R is the ratio of cross sections for p production 
by longitudinal and transverse photons. The nucleon helicities are omitted for brevity. 
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Element 


Measurement 


r 04 
'00 

Rer?4 

r 04 

'1-1 

r oo 

r ll 

Re r{ 

Im rf 

Im r\_ x 

/viO 

'oo 

Ax 

Rerf 

r\_ x 

Im rf 

Im rf_ x 


0.674 ±0.018 toml 
0.011 ±0.012 t° °ool 
-0.010 ± 0.013 ±g;gg| 
-0.058 ± 0.048 +°;°}? 
0.002 ± 0.034 ±°;°°<j 
-0.018 ±0.016 to \l 
0.122 ±0.018 +^oo5 
0.023 ± 0.016 +^oo9 
-0.119 ±0.018 Iggg 
0.093 ± 0.024 tamo 
0.008 ± 0.017 *j;°!g 
0.146 ±0.008 t° oot 
-0.004 ± 0.009 toMs 
-0.140 ±0.008 toiool 
0.002 ± 0.009 +^ooo 



Table 3: Spin density matrix elements for the elastic electroproduction of p mesons, measured 
for the 1996 data sample as the average values of the corresponding orthogonal functions of 
the p meson production and decay angles (see Appendix C of ref. [|30|]). The first errors are 
statistical, the second systematic. 
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Figure 6: Spin density matrix elements for elastic electroproduction of p mesons, measured 
for three values of Q 2 with the 1996 data sample. The inner error bars are statistical and the 
full error bars include the systematic errors added in quadrature. The dashed lines indicate the 
expected null values in the case of SCHC. 
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Figure 7: Spin density matrix elements for elastic electroproduction of p mesons, measured 
for three values of W with the 1996 data sample. The inner error bars are statistical and the 
full error bars include the systematic errors added in quadrature. The dashed lines indicate the 
expected null values in the case of SCHC. 
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Figure 8: Spin density matrix elements for elastic electroproduction of p mesons, measured for 
three values oft with the 1996 data sample. The inner error bars are statistical and the full error 
bars include the systematic errors added in quadrature. The dashed lines indicate the expected 
null values in the case of SCHC. 
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5.2 Helicity Conserving Amplitudes 



5.2.1 Ratio of the Longitudinal and Transverse Cross Sections 



After integration over the angles ip and (p, the angular distribution (eq. ^) takes the form 

W(cos6) oc 1 - r{g + (3 r$ - 1) cos 2 6 . 



(32) 



In Fig. 15, the cos# distributions for the 1996 data are presented for six bins in Q 2 , and the 



results of fits to eq. ( p2| ) are superimposed. As can be observed from the figures, the quality of 
the fits is good. The resulting measurements of r^ are in good agreement with those presented 
in Figs. |-| and in Tables |-§. 

In the case of SCHC, the matrix element Tqq provides a direct measurement of R, the ratio 
of cross sections for p production by longitudinal and transverse virtual photons (see Table |2|, 
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Figure 9: Acceptance corrected cos^ distributions for the 1996 data sample in six bins in Q 2 



The curves are the result of fits to the form of eq. (p2|). The errors on the data points are 
statistical only. 
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Figure 10: The ratio R of cross sections for elastic p meson electroproduction by longitudinal 
and transverse photons, measured in the SCHC approximation and presented as a function of 
Q 2 . For the present measurements (full circles), the inner error bars are statistical and the full 
error bars include the systematic errors added in quadrature. The other measurements are from 
HI [g§] and ZEUS Q in photoproduction, and from CHIO [g], NMC [§, E665 [§], HI ^ 
and ZEUS [0] in electroproduction. 
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(33) 



As the SCHC violating amplitudes are small compared to the helicity conserving ones (see 
section ^3j), eq. ( |33"1 ) can be used assuming SCHC to estimate R^\ 



The values of R deduced from eq. ( [I3"| ) using the results of the fits of the cos 6 distributions 
to eq. ( |32| ) are presented in Fig. [T0| (and in Table |7|) as a function of Q 2 , together with other 
measurements performed assuming SCHC [p]-|5|,[2T)|,[2"8|]. It is observed that R rises steeply at 
small Q 2 , and that the longitudinal cross section dominates over the transverse cross section 
f° r Q 2 ^ 2 GeV 2 . However, the rise is non-linear, with a weakening dependence at large Q 2 
values. No significant W dependence of the behaviour of R as a function of Q 2 is suggested by 
the comparison of the fixed target and HERA results. 

9 The Toi amplitude, which appears to be the dominant helicity-flip amplitude, corresponds to 8 ± 3% of the 
non-flip amplitudes ^/l^ool 2 + |2ii| 2 ( see section 5.3). A comparison of the forms of Tqq in columns 2 and 3 of 



Table indicates that the effect of SCHC violation on the measurement of R is 2.5 ± 1.5%. This is neglected. 
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5.2.2 Longitudinal— Transverse Interference 

In the case of NPEQ and SCHC, the decay angular distribution W(cos 9, if, <p) reduces to a 
function of two variables, cos# and ift, where ip = <p — Lp h the angle between the electron 
scattering plane and the p meson decay plane: 

W (cos 9,iff) = =- ( sin 2 9 (1 + e cos 2ip) 

87r 1 + e R I 

+ 2eR cos 2 9 - y/2e {1 + e) R cos 5 sin 29 cos ip\ . (34) 

Here 5 is the phase between the transverse T\\ and the longitudinal T 00 amplitudes: 

rooT^^lTooMTule-^ (35) 

and 

1 I , D 

cos 6 = -== (Re rf - Im r%) . (36) 

A two-dimensional plot of the cos 9 and ip variables is presented in Fig. |TT] for the 1996 



data. A fit of eq. (|34|) to these data gives: 



cos 5 = 0.925 ± 0.022 to oil ■ (3V) 



This number is in agreement within errors with the value of cos S computed from eqs. ( |33| ) and 
d3~6[ ) using the measurements of Tqq, Re rf and Im rf given in Table ^|. 

Fig. |T2](and Table ^ presents the measurements of cos 5 as a function of Q 2 , W and t. No 
significant evidence is found for a variation in the phase between the transverse and longitudinal 
amplitudes with these variables. That these amplitudes are nearly in phase was already observed 
at lower energy [|3l,|6|,[3~2"|]. 



5.2.3 The if) Distribution 

Fig. [T3| shows the distributions of the angle ip for five bins in Q 2 . They are well described by 
the function 

W(tl>) = — (1 + 2 e rj_! cos 2tfj) , (38) 

2n 

obtained from the integration over cos 9 of the function W(cos 9, if, <p) (eq. [Z3]), assuming 

SCHC. Measurements of the r\_ x matrix element extracted from fits to eq. ( |38| ) as a function 

of Q 2 , W and t are in good agreement with the measurements presented in Figs. ||— 1| and in 

Tables |]— 16|, which supports the fact that SCHC is a good approximation for the present data. 

10 The asymmetry P a between natural (<r N ) and unnatural {a u ) parity exchange can be determined, for trans- 
verse photons, from the measured matrix elements as: 

a N _ cr c/ 

P ° = n^ u = 0- + *R) (2r}_i - ^oo) , 
a + cr u 

and is found to be compatible with 1, as at lower energy [g,^]. This implies that NPE holds in the data at least 
for transverse photons. The measurement of the corresponding asymmetry for longitudinal photons would require 
two different values of e, i.e. two beam energies (see eq. (103) in ref. [|30|] >. 
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Figure 1 1 : Acceptance corrected plot of the event distribution in cos 9 and ip for the 1996 data 
sample. 
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Figure 12: Measurements of the cos 5 parameter as a function of Q 2 , W and t, obtained assum- 
ing SCHC and NPE from fits to the (cos 9, ip) distributions. The inner error bars are statistical 
and the full error bars include the systematic errors added in quadrature. 
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Figure 13: Acceptance corrected distributions of the ip angle for the 1996 data sample in five 
bins in Q 2 . The curves are the result of fits to the form of eq. ([38]). The errors on the data points 
are statistical only. 
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5.3 Helicity Flip Amplitudes 

Distributions of the angle are presented in Fig. [14] for six bins in Q 2 . These distributions, as 
well as the corresponding distributions for bins in W and t, exhibit significant variation in cos 0. 
Variation in cos 20 is compatible with zero. They are well described by the function 



W(<j>) oc 1 - £ cos 20 (2 rjj + r l m ) + yj2e{l + e) cos (2 r 5 n + r<> ) , (39) 

obtained from the integration of the decay angular distribution (eq. E31) over cos 9 and ip. 



The combinations of the matrix elements (2 r\ x + Tq ) and (2 rf x + r^), extracted as a 
function of Q 2 , VT and t, are presented in Fig. [15] (and in Table |9|). There is no indication for 
a significant deviation from zero of the combination (2 r\ x + r\ ), which is consistent with 
the measurements presented in Figs. |§— |] and in Tables 0—^]. In contrast, the combination 
(2 r\ x + Tq ) is significantly different from zero. As discussed in section pT] , this effect is 
attributed to a violation of SCHC for the matrix element Tq . 

As can be deduced from the second column of Table 0, the Tq matrix element is approx- 
imately proportional to the amplitude T i for a transverse photon to produce a longitudinal p 
meson: 

r °°- l + eR\T u \> (4U) 

where the term |T i| 2 has been neglected with respect to |T n | 2 in the denominator and the 



amplitudes T 00 and T 01 are assumed to be in phase and purely imaginary [31]]. 



With these approximations and with s ~ 1, the measurement of Tq allows the determination 
of the ratio of the T 01 amplitude to the non-flip amplitudes a/|T 00 | 2 + |T n | 2 for the present Q 2 
domain: 



l Toi l l Toi l ~ r s ./ 1 + ig (41) 

VlTooP + ITuP |Tn|Vr+fl °° V 2R 

~ 8 ±3%, (42) 

using the results in Table ^] and eq. (P3]). This value is of the order of magnitude, or slightly 
lower than those found, with large errors, at lower energy and for (Q 2 ) ~ 0.5 GeV 2 (15 — 20% 
for W ~ 2.5 GeV [§ and 11 - 14% for 10 < W < 16 GeV [§). 

The other helicity flip amplitudes are consistent with zero within the present measurement 
precision, as can be deduced from the fact that among the matrix elements which vanish under 
SCHC only the Tq element is measured to be non-zero. This is confirmed by the study of the if 
distribution. After integration over cos 9 and 0, the decay distribution (Q) reduces to 

W(ip) oc 1 - 2 rjij cos 2f . (43) 



31 




100 200 300 

* [deg.] 




100 200 300 

* [deg.] 




100 200 300 

* [deg.] 



© 

1 







60 


3.0 < Q 2 < 4.0 GeV^= 


40 


=^i^^f^ 


20 


-f^^r^ ^ 


n 


ii i i i i i i i i i i i i i ii 



100 200 300 

* [deg.] 



© 

T3 




100 200 300 

* [deg.] 



© 




100 200 300 

* [deg.] 



Figure 14: Acceptance corrected distributions of the <\> angle for the 1996 data sample in six 
bins in Q 2 . The curves are the result of fits to the form of eq. (|59"|). The errors on the data points 
are statistical only. 
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Figure 15: Measurements of the combinations of matrix elements 2r\ x + Tq and 2rf x + rQ , as 
a function of Q 2 , W and £, obtained from fits to the <p distributions. The inner error bars are 
statistical and the full error bars include the systematic errors added in quadrature. The dashed 
lines indicate the null values which are expected in the case of SCHC. 

This distribution is compatible with being constant for all bins in Q 2 , W and t, supporting 
the observation that the matrix element r°i x is consistent with zero. The expression for this 
matrix element contains a term proportional to T u Tf_ 1 , the interference between the helicity 
conserving transverse amplitude and the double-flip amplitude, and a term proportional to the 
square of the single flip contribution T x 2 (NPE is assumed). The constant (p distributions thus 
indicate that the helicity amplitudes Ti_i and T 10 are compatible with zero. 



,.0-i 
00 



Another way to study the amplitude Tio is to compare the measured values of the r\_ 1 and 
matrix elements, which are related by 



1-1 



> 



00 J 



(44) 



if, and only if, T 10 = (NPE is assumed). Relation (JHJ) is satisfied within errors for the 
measurements presented in Figs. |6|— 18| and in Tables |]— 1|. 



5.4 Comparison with Models 



Numerous models for the electroproduction of vector mesons based on VDM or QCD have 
been proposed. Most of them predict, for the present Q 2 domain, a linear increase with Q 2 
of the ratio R of the longitudinal to transverse cross sections, in disagreement with the results 
presented in Fig. [TO. 
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Figure 16: The ratio R of the longitudinal to transverse photon cross sections for elastic p meson 
electroproduction as a function of Q 2 . The data are the HERA measurements as in Fig. [K]. The 
curves are the predictions of the models of Royen and Cudell [ ]3~5| ] (solid), of Martin, Ryskin and 
Teubner (dashed) @ and of Schildknecht, Schuler and Surrow [f33|] (dotted), for the HERA 
energy range. 



However, several recent models predict a slower increase of R at high Q 2 [ |3Tl , [33l -|35l], which 
corresponds better to the trend in the data. One of them offers in addition full predictions for the 



spin density matrix elements [pT]]. In the rest of this section, we concentrate on the comparison 
of these model predictions with the present measurements. 



5.4.1 Generalised Vector Dominance 



A calculation based on the Generalised Vector Dominance Model (GVDM) has been performed 
by Schildknecht, Schuler and Surrow [|33|]. It takes into account a continuous mass spectrum 
of vector meson states, with destructive interferences between neighbouring states. This leads 
to a non-linear Q 2 dependence for the ratio R, in contrast with the conventional VDM predic- 
tions. The ratio R tends asymptotically to a constant value, defined by effective transverse and 
longitudinal masses which must be obtained from a fit to experimental data. The domain of 
applicability of the model extends in Q 2 down to photoproduction. 

In Fig. [Tj], the prediction of this model is compared to the measurement of R as a function of 
Q 2 , using the best set of parameters ("2-par. fit" in [|33|]). The data are the HERA measurements 
presented in Fig. [TO. 
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5.4.2 Parton-Hadron Duality 

Martin, Ryskin and Teubner have observed that QCD calculations of the p cross section that 
convolute the scattering amplitude with the p wave function give transverse cross sections which 
fall off too quickly with increasing Q 2 and thus lead to values of R which are too large at high 
Q 2 [J34]] • They have proposed an alternative approach, in which open qq production is considered 
in a broad mass interval containing the p meson. Hadronisation proceeds predominantly into 
two pion states, following phase space considerations. The hard interaction is modelled through 
two gluon exchange (or a gluon ladder), which induces a dependence on the parameterisation 
of the gluon density in the proton. The main uncertainties of the model come from the higher 
order corrections and from the choice of the mass interval embracing the p meson. However, 
the prediction for the ratio R of the cross sections has little sensitivity to these uncertainties. 

Fig. [R] presents the prediction of the model of Martin, Ryskin and Teubner, computed with 
the MRS(R4) parameterisation [f36|,[37l] for the gluon content of the proton. 



5.4.3 Quark Off-Shellness Model 

Another model based on lowest-order perturbative QCD calculations has been proposed by 
Royen and Cudell p5|]. The p meson production is computed from the qq Fock state of the 



photon, convoluted with the amplitude for hard scattering modelled as two-gluon exchange. A 
proton form factor and a meson vertex wave function, including Fermi motion, are part of the 
calculation. The specific feature of the model is that the constituent quarks are allowed to go 
off-shell. The W dependence of the cross section is not predicted, but the Q 2 and t dependences 
are. The uncertainties of the model come from the choice of the constituent quark mass m q and 
the Fermi momentum pp. 



The prediction of the model of Royen and Cudell is shown in Fig. |16|for m q = 0.3 GeV and 
Pf = 0.3 GeV. When m q and p F are varied by ± 50 MeV, the R value changes by about 15% 
and 30%, respectively, for Q 2 = 10 GeV 2 . 



5.4.4 Predictions of Polarisation 

Ivanov and Kirschner have provided predictions for the full set of 15 elements of the spin density 
matrix, based on perturbative QCD [plj]. This model predicts a violation of SCHC at high Q 2 , 
the largest helicity-flip amplitude being T i, with: 

|Too| > |2\i| > |T 01 | > |T 10 | > |Ti_i| (45) 

for the HERA kinematical domain. The ratios |Tn|/|T o|, |T i|/|T o| and |Ti |/|T o| depend on 
t, Q 2 , M and 7, where M is the invariant mass of the qq pair and 7 is the anomalous dimension 
of the gluon density (xg(x,Q 2 ) oc Q 2 " 1 ). The ratio |T 1 _ 1 |/|T 00 | depends also on the gluon 
density at the scale Q 2 /4. 
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Figure 17: Spin density matrix elements for elastic electroproduction of p mesons, for three 
values of Q 2 . The data are the same as in Fig. |5|. The curves are the predictions of the model 
of Ivanov and Kirschner [ ]3~T| ] for the GRV 94HO parameterisation of the gluon density in the 
proton. 
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Fig. Wj\ shows the predicted values of the matrix elements obtained with the parameterisa- 



tion GRV 94HO of the gluon density in the proton Q3q, p8|], compared to the measurements 
presented in Fig. |6|. This density is assumed to be valid throughout the range of Q 2 of the data. 
For higher Q 2 values, other parameterisations give predictions differing by much less than the 
measurement uncertainties. Reasonable agreement of the model predictions with the data is ob- 
served, with a correct prediction of the hierarchy between the amplitudes which are measured 
to be non-zero, and of the magnitude of the matrix element Tq . 



6 Cross Sections 

6.1 t Dependence of the ep Cross Section 

The acceptance corrected t distributions of the selected events with \t\ < 0.5 GeV 2 are presented 
in Fig. [T8| for five bins in Q 2 . To study the t dependence of elastic p production, these distri- 



butions are fitted as the sum of three exponentials corresponding to the elastic component, the 
diffractive component with proton dissociation and the non-resonant two-pion background. The 
elastic component is fitted with a free slope parameter b, whereas the contribution of diffractive 
p events with proton dissociation, which amounts to 11 ± 5% of the elastic signal, has a fixed 
slope parameter b p d — 2.5 ± 1.0 GeV~ 2 (see section p.2.2[ ).P1 The non-resonant background, 
amounting to 1 ± 1% of the signal, also has a fixed slope parameter, b nr = 0.3 ±0.1 GeV -2 , 
extracted from the present data at large \t\ values. 

The fitted exponential slope parameters, 6, for elastic events are presented as a function of 
Q 2 in Fig. [T9| (and in Table $T$>, together with HI [JT|,^§, ZEUS [g[Z^|3]5p and fixed target ^§\ 
measurements. The systematic errors are computed by varying the parameters of the Monte 
Carlo simulation used for the acceptance corrections (see section |3~TT| ), by varying the amounts of 
background contributions and their slopes within the quoted errors, and by varying the binning 
and the limits of the fits. 

The present measurements confirm the decrease of b when Q 2 increases from photoproduc- 
tion to the deep-inelastic domain, presumably reflecting the decrease of the transverse size of 
the virtual photon. It is also observed in Fig. [19| that at low Q 2 (Q 2 < 2 GeV 2 ) measurements at 
HERA lie systematically above the low energy fixed target results. This may indicate shrinkage 
of the diffractive peak as W increases. At higher Q 2 , given the experimental errors, no signif- 
icant information on a possible shrinkage of the t distribution can be extracted within the W 
range of the present experiment. 

The Q 2 evolution of the t distribution in the model of Royen and Cudell [551] is compared in 



Fig. [Z0|to the present measurements, in the form of the variable l/(|t|), which coincides with b 
for an exponential distribution. Q The trend of the data is reproduced. 

1 ' It should be noted that the slope parameter for low mass excited proton states could be larger than in the 
high mass region, from which the parameter b p d is extracted. The corresponding uncertainty is covered by the 
systematic errors quoted below. 

12 For the ZEUS measurements, the definitions of the slope differ somewhat: in the photoproduction case |[28||. 
the exponent of the t distribution was parameterised in a parabolic form, and only the linear term is plotted here; 
the fit in iQ was restricted to |*| < 0.4 GeV 2 and that in [g] was performed for |i| < 0.3 GeV 2 . 

13 In the present kinematic domain, the integration limits of \t\, \t\ min and |i| maa; , are such that \t\ min fv <C 
(|t|)and|i| m<KC >(|t|). 
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Figure 18: Acceptance corrected t distributions for the 1996 data sample, for five bins in 
Q 2 . The full curves correspond to a fit of the distributions as the sum of three exponentials, 
corresponding to the elastic signal (dashed curves), 11% background of proton dissociation 
events with slope b p d = 2.5 GeV -2 (dotted), and 1% non-resonant background with slope b nr = 
0.3 GeV 2 (not visible on the plots). The errors on the data points are statistical only. 
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Figure 19: Measurement of the slope parameter b of the exponential t dependence for elastic 
p production. For the present measurements (full circles), the inner error bars are statistical 
and the full error bars include the systematic errors added in quadrature. The other measure- 
ments are from HI |JZ§] and ZEUS [gS^ in photoproduction, and from CHIO [gfl, NMC 0], 
E665 |Q], HI [[]]] and ZEUS 0] in electroproduction. It should be noted that the definition of 
the parameter b is not unique (see text). 
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Figure 20: Q 2 dependence of the slope parameter b for p elastic production by HI (these and 
previous measurements [jl],^6|]), compared to the predictions of the model of Royen and Cud- 



ell [ p5Q for the HERA energy range, presented in the form of the variable l/(|t|). 
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6.2 Q 2 Dependence of the 7*p Cross Section 

The 7*p cross section for p elastic production is extracted from the ep cross section using the 
relation: 

d2(Ti dydQi PP) = r a(7> "* PV) = r aAl * P ~* PP) (1 + ei2) ' (46) 

where T is the flux of virtual photons [fTO|], given by: 

„ a(l-?/ + ?/ 2 /2) 



vyQ 2 



(47) 



a being the fine structure constant. The flux is integrated over each kinematic domain using the 
measured Q 2 and W dependences of the YP cross section. 

The 7*p cross section is presented in Fig. ^T] (and in Table [TT]) as a function of Q 2 , for a 
common value W = 75 GeV. It is obtained from the fits described in section ^, which take into 
account the Q 2 dependent skewing of the m^ mass distribution. The cross section is quoted 
for a relativistic Breit-Wigner distribution of the p mass, described by eqs. ( |T9| ) and (|20[), for 
the mass interval 

2 m^ < "Vtt < vn p + 5 T p . (48) 

The use of two alternative forms to eq. ( |20| ) in parameterising the width r(m 7r7r ) [ p4| ] would 
cause an increase of the cross section by 5%, which is included in the systematic errors. The 
background contributions of p diffractive production with proton dissociation, of uj and <p elastic 
production, and the non-resonant background are subtracted assuming the same distribution in 
Q 2 as for the signal. The uncertainties in these backgrounds are included in the systematic 
errors. The Q 2 dependent losses induced by the \t\ < 0.5 GeV 2 cut are corrected for on a bin-by- 
bin basis, according to the measured b slope parameters (see section |Q| ). The data are corrected 
for the losses of events due to noise in the detectors FMD and PRT (5 ± 3%) and LAr (10 ± 3%). 
Acceptance and efficiency effects and their errors are determined as described in section |3~1] . 
The errors on the extrapolations of the cross sections to the common value W = 75 GeV and 
to the quoted Q 2 values are estimated by varying the assumed W and Q 2 dependences of the 
cross section according to the limits of the present measurements. The radiative corrections are 
very small for the chosen value of the E — p z cut and for the procedure used to compute the 
kinematic variables (see section 1272]); an error of 4% accounts for the relevant uncertainties in 



the Q 2 and W dependences of the cross section, for higher order processes, and for detector 
effects not simulated in detail. The systematic errors on the cross section measurements also 
include an uncertainty of 2% in the luminosity, and the uncertainties due to limited Monte Carlo 
statistics. 

A parameterisation of the Q 2 dependence of the cross section in the form 
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Figure 21: Cross section measurements for the process YP ~^ PP plotted as a function of Q 2 
for W = 75 GeV (the ZEUS measurements [JZp have been scaled to W = 75 GeV. The inner 
error bars are statistical and the full error bars include the systematic errors added in quadrature. 
The curve corresponds to a fit to the present data of the form of eq. (|49|), with n = 2.24. 
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Figure 22: Cross section for the process 7*p — > pp, plotted as a function of Q 2 for W = 75 GeV 
(the data are the same as in Fig. ^1]). Photoproduction measurements by HI [[26p and ZEUS [ |28| ] 
are also shown. The curves are the predictions of the models of Royen and Cudell [ |35| ] (solid), 
of Martin, Ryskin and Teubner [ |3"4"| ] (dashed) and of Schildknecht, Schuler and Surrow p3| ] 
(dotted). 



is shown superimposed on Fig. [2_y. It is obtained by a fit to the present data with the result 

n = 2.24 ±0.09. (50) 

The uncertainty on this value is determined using the statistical and the non-correlated system- 
atic errors only. The nominal normalisations are used for the 1995 and 1996 data sets, which 
agree within one standard deviation. The quality of the fit for the full Q 2 range 1 < Q 2 < 
35 GeV 2 is good: rVridf =13.3/20. 



Fig. fTi] presents in addition the measurements of the ZEUS collaboration [f2p, scaled to the 
value W = 75 GeV. Agreement is observed between the results of the two experiments. 



In Fig. [22|, the Q 2 dependence of the YP cross section, including photoproduction measure- 
ments p5], [ZS1], is compared with the predictions of the models of Schildknecht, Schuler and 
Surrow [[331], °f Martin, Ryskin and Teubner [ |34| ] and of Royen and Cudell [p3|]. The latter 
model describes the data well down to the photoproduction region. 
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6.3 W Dependence of the 7*p Cross Section 

The 7*p cross section for p elastic production is presented as a function of W for six values of 
Q 2 in Fig. ^3] (and in Table |12|). The extrapolations of the measured cross sections to the chosen 
Q 2 values are performed using the Q 2 dependence given by eqs. ( ft9| ) and (|5TJ[). Corrections and 
systematic errors are determined as described in section |6\2| . 



To quantify the W dependence of the cross section, a fit is performed for each Q 2 bin to a 
power law: 

a(i*p)<xW s , (51) 

as shown in Fig. ^3|. Only the statistical and the non-correlated systematic errors are used in the 
fits, and the values of x 2 /ndf are reasonable for all Q 2 bins. 

In a Regge context, the parameter 5 can be related to the exchange trajectory :Q 

S ~ 4 [a((t)) - 1] . (52) 

The trajectory is assumed to take a linear form: 

a{t) = a(0) + at. (53) 

To extract the effective trajectory intercept a(0), (\t\) = lib is taken from the measured 
values (see section |67fl). In the absence of a measurement of the Q 2 dependence of the shrinkage 



of the t distribution with increasing W, the value a' == 0.25 GeV -2 is assumed, as measured 
in hadron— hadron interactions [|4l|]. The values obtained for the intercept a(0) as a function of 
Q 2 are shown in Fig. 123] (and in Table O). The inner error bars come from the statistical and 



non-correlated systematic uncertainties on the cross section measurements. The sensitivity to 
the choice of a' is shown by the outer bars, which contain the variation due to the assumption 
a' = (i.e. no shrinkage) added in quadrature. The measurements are compared to the values 
1.08— 1 . 10 obtained from fits to the total and elastic hadron-hadron cross sections pT| , ^2| ] . They 
suggest that the intercept of the effective trajectory governing high Q 2 p electroproduction is 
larger than that describing elastic and total hadronic cross sections. 

It should be noted that several studies (see e.g. [p3|]) indicate that QCD based predictions 
for the W dependence of the 7*p cross section are affected by large uncertainties. These are 
related particularly to the assumptions made concerning the appropriate factorisation scale and 
the p wave function, and also to the choice of the parameterisation of the gluon distribution in 
the proton. 



14 Strictly speaking, this applies if the W dependence of the integrated cross section J da /At At is the same, 
over the relevant W domain, as the W dependence of the differential cross section Au/At for t = (t). 

43 



■§ 



Q. 

t 
P. 



10 



10 



10 




40 50 60 70 80 90100 



W [GeV] 



Figure 23: Cross section for the process Y'p ~~ * PP as a function of W 7 for several values of Q 2 . 
The inner error bars are statistical and the full error bars include the systematic errors added in 
quadrature. The lines correspond to a fit of the form of eq. (|5"T|). 
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Figure 24: Q 2 dependence of the intercept a(0) (see eqs. |5T] — |53|). The inner error bars rep 



resent the statistical and non-correlated systematic uncertainties on the cross section measure- 
ments, the outer error bars include the variation of the intercept a(0) when assuming a' = 0, 
added in quadrature. The dashed lines represent the range of values obtained for the "soft 
pomeron" intercept, as derived from fits to total and elastic hadron— hadron cross section mea- 



surements [41, 42] 



7 Summary and Conclusions 

The elastic electroproduction of p mesons has been studied at HERA with the HI detector, for 
1 < Q 2 < 60 GeV 2 and 30 < W < 140 GeV. 

The shape of the (nn) mass distribution has been studied as a function of Q 2 . It indicates 
significant skewing at low Q 2 , which gets smaller with increasing Q 2 . 

The full set of 15 elements of the p spin density matrix has been measured as a function 
of Q 2 , W and t, using the decay angular distributions defined in the helicity frame. Except for 
a small but significant deviation from zero of the Tq element, s-channel helicity conservation 
is found to be a good approximation. For Q 2 > 2 GeV 2 , the longitudinal Y'p cross section 
becomes larger than the transverse cross section, and the ratio R reaches the value R ~ 3 for Q 2 
~ 20 GeV 2 . The phase 5 between the longitudinal and transverse amplitudes is measured to be 
cos 5 = 0.93 ± 0.03, assuming natural parity exchange and s-channel helicity conservation. The 
dominant helicity flip amplitude T A =o,a =i is found to be 8 ± 3% of the non-flip amplitudes. 
A model based on GVDM [ J3~3"l ] and models based on perturbative QCD [ |34| , |35l ] reproduce the 
flattening of the ratio R observed at high Q 2 . A QCD based prediction pl| ] is in qualitative 
agreement with the measurement of the 15 matrix elements, in that it reproduces the observed 
hierarchy between the amplitudes which are measured to be non-zero and the magnitude of the 
matrix element Tq . 

The t distribution for p electroproduction has been studied and the exponential slope pa- 
rameter b is found to decrease when Q 2 increases from photoproduction to the deep-inelastic 
domain. 
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The 7*p — » pp cross section has been measured over the domain 1 < Q 2 < 35 GeV 2 and 
follows a Q 2 dependence of the form l/(Q 2 + rn 2 ) n , with n = 2.24 ± 0.09. This dependence 
is well described by a model based on QCD [p~~ 



The W dependence of the 7*p — > pp cross section has been measured for six values of Q 2 . 
The measurements suggest that the intercept of the effective trajectory governing high Q 2 p 
electroproduction is larger than that describing elastic and total hadronic cross sections. 
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Element 



2.5 < Q 2 < 3.5 GeV 2 



3.5 < Q 2 < 6.0 GeV" 



6.0 < Q' 2 < 60 GeV 2 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 



r 04 
'00 

Rer^ 

r 04 
' 1-1 



00 



'11 

Re r\ 



'i-i 
Im r 2 

Im r^_ 1 



oo 



'ii 
Rerf 



'i-i 
Imrfo 

Im r^_ x 



0.639 ± 

0.018 ± 

-0.020 ± 

-0.011 ± 

-0.019 ± 

0.003 ± 

0.147 ± 

0.006 ± 

-0.156 ± 

0.099 ± 

0.002 ± 

0.149 ± 

-0.019 ± 

-0.124 ± 

0.022 ± 



0.031 
0.020 
0.023 
0.081 
0.057 
0.028 
0.032 
0.028 
0.032 
0.040 
0.029 
0.013 
0.016 
0.013 
0.016 



+0.013 
-0.010 
+0.004 
-0.004 
+0.002 
-0.003 
+0.013 
-0.022 
+0.016 
-0.003 
+0.016 
-0.003 
+0.013 
-0.007 
+0.008 
-0.004 
+0.014 
-0.010 
+0.016 
-0.002 
+0.002 
-0.007 
+0.002 
-0.002 
+0.005 
-0.001 
+0.001 
-0.001 
+0.001 
-0.002 



0.695 ± 

-0.019 ± 

-0.020 ± 

-0.085 ± 

0.021 ± 

-0.022 ± 

0.103 ± 

0.049 ± 

-0.098 ± 

0.081 ± 

0.005 ± 

0.142 ± 

0.004 ± 

-0.146 ± 

-0.014 ± 



0.031 
0.020 
0.022 
0.082 
0.057 
0.028 
0.031 
0.028 
0.031 
0.041 
0.029 
0.013 
0.016 
0.013 
0.016 



+0.019 
-0.018 
+0.009 
-0.003 
+0.008 
-0.001 
+0.021 
-0.013 
+0.001 
-0.006 
+0.015 
-0.016 
+0.003 
-0.005 
+0.015 
-0.015 
+0.018 
-0.007 
+0.009 
-0.013 
+0.006 
-0.006 
+0.003 
-0.002 
+0.006 
-0.006 
+0.004 
-0.005 
+0.003 
-0.004 



0.748 ± 0.033 toml 
0.036 ± 0.022 +°;°^ 
0.016 ± 0.023 +°;°^ 

-0.078 ± 0.092 1°;°^ 
0.010 ± 0.063 +£009 

-0.042 ± 0.030 tomr 

-0.010 
-0.013 



0.081 ±0.031 



0.007 ± 0.030 tofA 
-0.067 ± 0.032 t °oll 

0.107 ± 0.047 1^04 
0.014 ± 0.032 tao?2 

0.133 ± 0.014 ;°;°04 
0.003 ± 0.016 +H?o 

+0.004 
-0.002 
+0.006 
-0.003 



-0.146 ±0.014 
-0.001 ± 0.016 



Table 4: Spin density matrix elements for elastic electroproduction of p mesons, measured 
for three values of Q 2 with the 1996 data sample. The first errors are statistical, the second 
systematic. 
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Element 



40 < W < 60 GeV 



60 < W < 80 GeV 



80 < W < 100 GeV 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 



r, 



01 
00 



Rer°4 

r 04 
' 1-1 



00 



' 11 
Re r\ 



'i-i 

Im r? 

Im rf^ 



oo 



' n 
Rerfo 



' i-i 
Im r^ 

Im ri_ x 



0.671 ±0. 
-0.011 ±0. 
-0.021 ± 0. 
-0.048 ± 0. 
-0.013 ± 0. 
-0.002 ± 0. 

0.225 ± 0. 
-0.030 ± 0. 
-0.132 ±0. 

0.030 ± 0. 
-0.009 ± 0. 

0.177 ±0. 
-0.014 ± 0. 
-0.148 ±0. 
-0.005 ± 0. 



+0.031 
-0.025 
+0.010 
-0.007 
+0.005 
-0.006 
+0.021 
-0.019 
+0.008 
-0.009 
+0.006 
-0.006 
+0.002 
-0.005 
+0.012 
-0.008 
+0.016 
-0.010 
+0.003 
-0.003 
nOQ +0.004 
VAy -0.001 
O1 o +0.015 
VLJ -0.009 
01fi +0.004 
U1U -0.008 
01 o +0.003 
V1J -0.004 
01fi +0.002 
U1U -0.002 



031 
020 
023 
081 
057 
028 
031 
028 
032 
041 



0.719 ± 
0.025 ± 
0.000 ± 

-0.151 ±0 
0.080 ± 

-0.018 ±0 
0.113 ±0 
0.105 ±0 

-0.151 ±0 
0.192 ±0 

-0.015 ± 
0.118 ±0 

-0.005 ± 

-0.160 ±0 
0.021 ± 



mi +0.051 
VJ1 -0.033 

020 +oon 

.VZ.V _ .009 
090 +0.011 
WZJ -0.010 

'•082 ±K 

nc 7 +0.003 

UJ ' -0.006 

. n oo +0.013 

'.VZ.0 _ .022 

1 031 +0 - 007 
• WJ1 -0.010 

. m o +0.006 

i.UZS _ 008 
, 03 1+0.013 

■ JJ1 -0.009 
i041 +0.033 
• W1 -0.013 

09 q +0.004 
.vz.y _ .oio 

on +0.011 
• uu -0.009 
01fi +0.007 
.U1U _ .002 

n ,o +0.005 
■ uu -0.006 

01fi +0.002 
.W1U _ .002 



0.687 ± 

0.052 ± 

-0.028 ± 

0.043 ± 

-0.060 ± 

-0.023 ± 

0.083 ± 

0.032 ± 

-0.068 ± 

0.114 ± 

0.027 ± 

0.125 ± 

-0.018 ± 

-0.115 ± 

-0.014 ± 



0.033 
0.021 
0.024 
0.089 
0.062 
0.030 
0.033 
0.030 
0.033 
0.045 
0.031 
0.014 
0.017 
0.014 
0.017 



+0.031 
-0.016 
+0.006 
-0.003 
+0.010 
-0.004 
+0.029 
-0.026 
+0.021 
-0.024 
+0.022 
-0.027 
+0.052 
-0.044 
+0.009 
-0.006 
+0.020 
-0.024 
+0.007 
-0.002 
+0.007 
-0.008 
+0.009 
-0.009 
+0.007 
-0.006 
+0.013 
-0.031 
+0.009 
-0.011 



Table 5: Spin density matrix elements for elastic electroproduction of p mesons, measured 
for three values of W with the 1996 data sample. The first errors are statistical, the second 
systematic. 



0.1 < Itl < 0.2 GeV" 



Element 



0.0 < Itl < 0.1 GeV 2 



0.2 < |t| < 0.5 GeV 2 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 



r 04 
'00 

Rer ^ 

r 04 
'1-1 



00 



'11 
Re r\ 



'1-1 
Im r 2 10 

Im r 2 _ 1 



00 



'11 
Rerfo 



'1-1 
Imrf 

Im r\_ x 



0.686 ± 

0.010 ± 

-0.011 ± 

-0.083 ± 

0.016 ± 

-0.032 ± 

0.098 ± 

0.020 ± 

-0.136 ± 

0.090 ± 

-0.003 ± 

0.155 ± 

-0.021 ± 

-0.143 ± 

0.004 ± 



0.031 
0.020 
0.023 
0.083 
0.058 
0.028 
0.031 
0.028 
0.031 
0.041 
0.029 
0.013 
0.016 
0.013 
0.016 



+0.050 
-0.039 
+0.009 
-0.014 
+0.012 
-0.005 
+0.026 
-0.031 
+0.018 
-0.016 
+0.022 
-0.013 
+0.013 
-0.015 
+0.014 
-0.010 
+0.007 
-0.003 
+0.055 
-0.038 
+0.015 
-0.027 
+0.005 
-0.011 
+0.014 
-0.007 
+0.005 
-0.006 
+0.002 
-0.000 



0.706 ± 

0.021 ± 

-0.011 ± 

-0.005 ± 

0.003 ± 

-0.044 ± 

0.134 ± 

0.045 ± 

-0.143 ± 

0.069 ± 

0.015 ± 

0.138 ± 

0.014 ± 

-0.122 ± 

-0.002 ± 



0.031 
0.020 
0.022 
0.082 
0.057 
0.028 
0.030 
0.028 
0.031 
0.041 
0.029 
0.013 
0.016 
0.013 
0.016 



+0.059 
-0.041 
+0.009 
-0.002 
+0.013 
-0.009 
+0.021 
-0.008 
+0.007 
-0.007 
+0.007 
-0.018 
+0.006 
-0.012 
+0.005 
-0.008 
+0.011 
-0.007 
+0.033 
-0.012 
+0.007 
-0.012 
+0.009 
-0.010 
+0.001 
-0.007 
+0.004 
-0.006 
+0.005 
-0.003 



0.634 ± 

-0.001 ± 

-0.005 ± 

-0.058 ± 

-0.030 ± 

0.029 ± 

0.170 ± 

0.023 ± 

-0.078 ± 

0.132 ± 

0.012 ± 

0.138 ± 

0.003 ± 

-0.152 ± 

0.001 ± 



0.033 
0.021 
0.024 
0.087 
0.061 
0.030 
0.033 
0.031 
0.033 
0.044 
0.031 
0.014 
0.017 
0.014 
0.017 



+0.042 
-0.032 
+0.013 
-0.005 
+0.005 
-0.012 
+0.015 
-0.006 
+0.006 
-0.006 
+0.010 
-0.013 
+0.014 
-0.009 
+0.009 
-0.005 
+0.025 
-0.007 
+0.020 
-0.050 
+0.012 
-0.013 
+0.012 
-0.002 
+0.003 
-0.007 
+0.004 
-0.001 
+0.004 
-0.004 



Table 6: Spin density matrix elements for elastic electroproduction of p mesons, measured for 
three values oft with the 1996 data sample. The first errors are statistical, the second systematic. 
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Q 2 (GeV 2 ) 


R 


= <JlI 


7 T 


1.8 


1.03 


+0.16 
-0.16 


+0.10 
-0.10 


2.7 


1.75 


+0.36 
-0.29 


+0.30 
-0.28 


3.4 


2.25 


+0.42 
-0.34 


+0.13 
-0.10 


4.8 


2.22 


+0.46 
-0.36 


+0.14 
-0.06 


7.2 


2.67 


+0.70 
-0.50 


+0.32 
-0.10 


10.9 


3.38 


+1.30 
-0.82 


+0.28 
-0.46 


19.7 


2.60 


+ 1.19 
-0.72 


+0.27 
-0.08 



Table 7: Measurement of the ratio R = ol/ctt for seven values of Q 2 , with (W) = 75 GeV, 
obtained from the measurement of the matrix element ?-qq, assuming SCHC. The first errors are 
statistical, the second systematic. 



Q 2 (GeV 2 ) 


W (GeV) 


\t\ (GeV 2 ) 


cos 5 


2.5 - 3.5 
3.5-6.0 
6.0 - 60. 


30 - 100 
30 - 120 
30 - 140 


0.0-0.5 
0.0-0.5 
0.0-0.5 


0.867 ± 0.051 tao?9 
0.841 ± 0.056 t H 
0.964 d= 0.071 ±g;gil 


2.5 - 60. 
2.5 - 60. 
2.5 - 60. 


40-60 
60-80 
80 - 100 


0.0-0.5 
0.0-0.5 
0.0-0.5 


0.922 ± 0.053 lg;{£° 
0.903 ± 0.064 +° f 2 2 8 
0.690 ± 0.101 tof 4 l 


2.5 - 60. 
2.5 - 60. 
2.5 - 60. 


30 - 140 

30 - 140 
30 - 140 


0.0-0.1 
0.1-0.2 
0.2-0.5 


o.9i5 d= 0.039 ±o:8ii 

0.904 ± 0.063 ±g:Sss 
0.868 ± 0.060 18;°jjjj 



Table 8: Measurements of the cos 5 parameter as a function of Q 2 , W and t, obtained under 
NPE and the SCHC approximation from fits to the (cos 9, ip) distributions. The first errors are 
statistical, the second systematic. 
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Q 2 (GeV 2 ) 


W (GeV) 


|*| (GeV") 


2 r}i+ r4 


2 r 5 n + rl 


2.5-3.0 
3.0-4.0 
4.0-6.0 
6.0-9.0 
9.0-14. 
14. - 60. 


30 - 100 
30 - 100 
30 - 120 
30 - 140 
30 - 140 
30 - 140 


0.0 - 0.5 
0.0 - 0.5 
0.0 - 0.5 
0.0 - 0.5 
0.0 - 0.5 
0.0 - 0.5 


0.046 ± 0.083 t f 5 9 

-0.140 ± 0.065 +°;°36 
-0.079 ± 0.072 toml 
-0.023 ± 0.084 to%il 
0.006 ±0.119+^ 
-0.173 ± 0.156 ±S:Jjg 


0.097 ± 0.039 toml 
0.115 ± 0.034 toml 
0.120 ± 0.036 lg;^ 
0.109 d= 0.043 ±8;8Si 

0.216 ±0.054 IS 
0.113 ± 0.077 ±g:o4o 


2.5 - 60.0 
2.5 - 60.0 
2.5 - 60.0 


40-60 
60-80 
80 - 100 


0.0 - 0.5 
0.0 - 0.5 
0.0 - 0.5 


-0.118 ± 0.066 lg;^ 
-0.040 ± 0.069 !{j;{£jj 
-0.106 ± 0.074 +oX 


0.025 ± 0.033 tomt 
0.175 ± 0.034 1{X 
0.183 ± 0.039 t°oml 


2.5 - 60.0 
2.5 - 60.0 
2.5 - 60.0 
2.5 - 60.0 


30 - 140 
30 - 140 
30 - 140 
30 - 140 


0.0-0.1 
0.1-0.2 
0.2 - 0.3 
0.3 - 0.5 


-0.060 ± 0.049 t ml 
0.012 ± 0.068 ±8:055 

-0.053 ± 0.090 +° ° H 
-0.182 ± 0.085 ±8;8Ii 


0.092 ± 0.025 +° ° l 8 
0.114 ± 0.033 +g;gj| 
0.126 ± 0.044 ±°;°£ 
0.196 ± 0.046 !jj:jgg 



Table 9: Measurements of the combinations of matrix elements 2rj 



ii 



?o and 2rj > 1 



±r, 



oo 



as a 



function of Q 2 , W and t, obtained from fits to the <fi distributions. The first errors are statistical, 
the second systematic. 



Q 2 (GeV 2 ) 


b (GeV-') 


1.8 


8.0 


±0.5 ±8:1 


3.1 


7.1 


±0.4 t° i 


4.8 


5.5 


±0.5 t° i 


7.2 


6.2 


±0.6 I** 


10.9 


5.6 


±0.8 ±ja 


19.7 


4.7 


±1.0 18:? 



Table 10: Measurement of the slope parameter b of the exponential t dependence for six values 
of Q 2 , with (W) = 75 GeV. The first errors are statistical, the second systematic. 
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Q 2 (GeV 2 ) 


a(Y 


p -> pp) (nb) 


1.1 


2129 


±369 


+275 
-275 


1.4 


1610 


± 194 


+207 
-207 


1.7 


1186 


±155 


+153 
-153 


2.3 


681 


±83 


+88 
-88 


2.7 


432 


±39 


+46 
-31 


3.0 


399 


±34 


+42 
-31 


3.3 


314 


±29 


+41 
-32 


3.8 


261 


±24 


+27 
-24 


4.2 


206 


±21 


+23 

-24 


4.7 


157 


± 17 


+14 
-17 


5.3 


120 


± 14 


+14 
-17 


6.0 


106 


±13 


+9 
-9 


6.7 


79 


±10 


+11 

-7 


7.5 


81 


±10 


+7 
-11 


8.4 


50.7 


±7.3 


+4.4 
-3.7 


9.4 


47.5 


±6.7 


+4.0 
-4.9 


10.9 


27.5 


±4.1 


+2.5 
-2.5 


13.0 


19.9 


±3.1 


+1.6 
-1.6 


15.4 


17.7 


±3.3 


+1.9 
-1.7 


18.3 


11.6 


±2.7 


+1.3 
-1.3 


22.8 


6.0 


±1.5 


+0.8 
-0.6 


35.0 


1.6 


±0.5 


+0.2 
-0.2 



Table 11: Measurement of the cross section for the process ^*p — > pp as a function of Q 2 for 
W = 75 GeV. The first errors are statistical, the second systematic. 
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Q 2 {GeV z ) 


W (GeV) 


a(r 


'p — ► pp) (nb) 


2.0 


49 


718 


±85 


+92 
-92 




65 


991 


± 118 


+128 
-128 




86 


1025 


± 117 


+133 
-133 




116 


1002 


± 118 


+129 
-129 


3.1 


40 


296 


±24 


+27 
-20 




60 


318 


±28 


+27 
-21 




80 


410 


±34 


+34 
-28 


4.8 


40 


125 


±13 


+11 
-11 




60 


137 


±16 


+ 12 
-10 




80 


160 


±19 


+16 
-15 




100 


168 


±21 


+15 
-14 


7.2 


50 


60.3 


±7.9 


+5.5 
-3.9 




70 


76.2 


±10.8 


+8.9 
-5.7 




90 


96.7 


± 14.3 


+12.2 
-6.2 




110 


94.7 


±15.0 


+8.4 
-6.1 




130 


76.9 


±18.8 


+11.8 
-6.8 


10.9 


50 


17.9 


±3.5 


+1.7 
-1.9 




70 


34.6 


±5.7 


+3.7 
-2.2 




90 


38.9 


±7.1 


+3.0 
-3.2 




110 


30.1 


±6.8 


+2.4 
-2.0 




130 


42.5 


±9.7 


+3.5 
-2.4 


19.7 


60 


7.2 


±1.6 


+0.6 
-0.6 




80 


10.4 


±2.3 


+0.7 
-1.0 




100 


9.4 


±2.4 


+1.0 
-0.7 




120 


14.9 


±3.2 


+1.1 
-1.1 



Table 12: Measurement of the cross section for the process 7*p — »■ pp as a function of W, for 
several values of Q 2 . The first errors are statistical, the second systematic. 
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Q 2 (GeV 2 ) 


a(0) 


2.0 


1.13 


±0.05 


+0.00 
-0.03 


3.1 


1.15 


±0.04 


+0.00 
-0.04 


4.8 


1.12 


±0.04 


+0.00 
-0.04 


7.2 


1.15 


±0.06 


+0.00 
-0.04 


10.9 


1.23 


±0.06 


+0.00 
-0.05 


19.7 


1.27 


±0.11 


+0.00 
-0.05 



Table 13: Measurements of the a(0) parameter (see eqs. I5T] — p31) as a function of Q 2 . The first 
error represents the statistical and non-correlated systematic uncertainties on the cross section 
measurements; the second error represents the variation due to the assumption a' = (i.e. no 
shrinkage). 
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